home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_rotg.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  48.3 KB  |  1,678 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_rotg (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    float a = -1.5f;
  13.    float b = -1.5f;
  14.    float c;
  15.    float s;
  16.    float r_expected = -2.12132034356f;
  17.    float z_expected = 1.41421356237f;
  18.    float c_expected = 0.707106781187f;
  19.    float s_expected = 0.707106781187f;
  20.    cblas_srotg(&a, &b, &c, &s);
  21.    gsl_test_rel(a, r_expected, flteps, "srotg(case 166)");
  22.    gsl_test_rel(b, z_expected, flteps, "srotg(case 167)");
  23.    gsl_test_rel(c, c_expected, flteps, "srotg(case 168)");
  24.    gsl_test_rel(s, s_expected, flteps, "srotg(case 169)");
  25.   };
  26.  
  27.  
  28.   {
  29.    float a = -1.5f;
  30.    float b = -1.0f;
  31.    float c;
  32.    float s;
  33.    float r_expected = -1.80277563773f;
  34.    float z_expected = 0.554700196225f;
  35.    float c_expected = 0.832050294338f;
  36.    float s_expected = 0.554700196225f;
  37.    cblas_srotg(&a, &b, &c, &s);
  38.    gsl_test_rel(a, r_expected, flteps, "srotg(case 170)");
  39.    gsl_test_rel(b, z_expected, flteps, "srotg(case 171)");
  40.    gsl_test_rel(c, c_expected, flteps, "srotg(case 172)");
  41.    gsl_test_rel(s, s_expected, flteps, "srotg(case 173)");
  42.   };
  43.  
  44.  
  45.   {
  46.    float a = -1.5f;
  47.    float b = -0.1f;
  48.    float c;
  49.    float s;
  50.    float r_expected = -1.50332963784f;
  51.    float z_expected = 0.0665190105238f;
  52.    float c_expected = 0.997785157857f;
  53.    float s_expected = 0.0665190105238f;
  54.    cblas_srotg(&a, &b, &c, &s);
  55.    gsl_test_rel(a, r_expected, flteps, "srotg(case 174)");
  56.    gsl_test_rel(b, z_expected, flteps, "srotg(case 175)");
  57.    gsl_test_rel(c, c_expected, flteps, "srotg(case 176)");
  58.    gsl_test_rel(s, s_expected, flteps, "srotg(case 177)");
  59.   };
  60.  
  61.  
  62.   {
  63.    float a = -1.5f;
  64.    float b = 0.0f;
  65.    float c;
  66.    float s;
  67.    float r_expected = -1.5f;
  68.    float z_expected = -0.0f;
  69.    float c_expected = 1.0f;
  70.    float s_expected = -0.0f;
  71.    cblas_srotg(&a, &b, &c, &s);
  72.    gsl_test_rel(a, r_expected, flteps, "srotg(case 178)");
  73.    gsl_test_rel(b, z_expected, flteps, "srotg(case 179)");
  74.    gsl_test_rel(c, c_expected, flteps, "srotg(case 180)");
  75.    gsl_test_rel(s, s_expected, flteps, "srotg(case 181)");
  76.   };
  77.  
  78.  
  79.   {
  80.    float a = -1.5f;
  81.    float b = 0.1f;
  82.    float c;
  83.    float s;
  84.    float r_expected = -1.50332963784f;
  85.    float z_expected = -0.0665190105238f;
  86.    float c_expected = 0.997785157857f;
  87.    float s_expected = -0.0665190105238f;
  88.    cblas_srotg(&a, &b, &c, &s);
  89.    gsl_test_rel(a, r_expected, flteps, "srotg(case 182)");
  90.    gsl_test_rel(b, z_expected, flteps, "srotg(case 183)");
  91.    gsl_test_rel(c, c_expected, flteps, "srotg(case 184)");
  92.    gsl_test_rel(s, s_expected, flteps, "srotg(case 185)");
  93.   };
  94.  
  95.  
  96.   {
  97.    float a = -1.5f;
  98.    float b = 1.0f;
  99.    float c;
  100.    float s;
  101.    float r_expected = -1.80277563773f;
  102.    float z_expected = -0.554700196225f;
  103.    float c_expected = 0.832050294338f;
  104.    float s_expected = -0.554700196225f;
  105.    cblas_srotg(&a, &b, &c, &s);
  106.    gsl_test_rel(a, r_expected, flteps, "srotg(case 186)");
  107.    gsl_test_rel(b, z_expected, flteps, "srotg(case 187)");
  108.    gsl_test_rel(c, c_expected, flteps, "srotg(case 188)");
  109.    gsl_test_rel(s, s_expected, flteps, "srotg(case 189)");
  110.   };
  111.  
  112.  
  113.   {
  114.    float a = -1.5f;
  115.    float b = 1.5f;
  116.    float c;
  117.    float s;
  118.    float r_expected = 2.12132034356f;
  119.    float z_expected = -1.41421356237f;
  120.    float c_expected = -0.707106781187f;
  121.    float s_expected = 0.707106781187f;
  122.    cblas_srotg(&a, &b, &c, &s);
  123.    gsl_test_rel(a, r_expected, flteps, "srotg(case 190)");
  124.    gsl_test_rel(b, z_expected, flteps, "srotg(case 191)");
  125.    gsl_test_rel(c, c_expected, flteps, "srotg(case 192)");
  126.    gsl_test_rel(s, s_expected, flteps, "srotg(case 193)");
  127.   };
  128.  
  129.  
  130.   {
  131.    float a = -1.0f;
  132.    float b = -1.5f;
  133.    float c;
  134.    float s;
  135.    float r_expected = -1.80277563773f;
  136.    float z_expected = 1.80277563773f;
  137.    float c_expected = 0.554700196225f;
  138.    float s_expected = 0.832050294338f;
  139.    cblas_srotg(&a, &b, &c, &s);
  140.    gsl_test_rel(a, r_expected, flteps, "srotg(case 194)");
  141.    gsl_test_rel(b, z_expected, flteps, "srotg(case 195)");
  142.    gsl_test_rel(c, c_expected, flteps, "srotg(case 196)");
  143.    gsl_test_rel(s, s_expected, flteps, "srotg(case 197)");
  144.   };
  145.  
  146.  
  147.   {
  148.    float a = -1.0f;
  149.    float b = -1.0f;
  150.    float c;
  151.    float s;
  152.    float r_expected = -1.41421356237f;
  153.    float z_expected = 1.41421356237f;
  154.    float c_expected = 0.707106781187f;
  155.    float s_expected = 0.707106781187f;
  156.    cblas_srotg(&a, &b, &c, &s);
  157.    gsl_test_rel(a, r_expected, flteps, "srotg(case 198)");
  158.    gsl_test_rel(b, z_expected, flteps, "srotg(case 199)");
  159.    gsl_test_rel(c, c_expected, flteps, "srotg(case 200)");
  160.    gsl_test_rel(s, s_expected, flteps, "srotg(case 201)");
  161.   };
  162.  
  163.  
  164.   {
  165.    float a = -1.0f;
  166.    float b = -0.1f;
  167.    float c;
  168.    float s;
  169.    float r_expected = -1.00498756211f;
  170.    float z_expected = 0.099503719021f;
  171.    float c_expected = 0.99503719021f;
  172.    float s_expected = 0.099503719021f;
  173.    cblas_srotg(&a, &b, &c, &s);
  174.    gsl_test_rel(a, r_expected, flteps, "srotg(case 202)");
  175.    gsl_test_rel(b, z_expected, flteps, "srotg(case 203)");
  176.    gsl_test_rel(c, c_expected, flteps, "srotg(case 204)");
  177.    gsl_test_rel(s, s_expected, flteps, "srotg(case 205)");
  178.   };
  179.  
  180.  
  181.   {
  182.    float a = -1.0f;
  183.    float b = 0.0f;
  184.    float c;
  185.    float s;
  186.    float r_expected = -1.0f;
  187.    float z_expected = -0.0f;
  188.    float c_expected = 1.0f;
  189.    float s_expected = -0.0f;
  190.    cblas_srotg(&a, &b, &c, &s);
  191.    gsl_test_rel(a, r_expected, flteps, "srotg(case 206)");
  192.    gsl_test_rel(b, z_expected, flteps, "srotg(case 207)");
  193.    gsl_test_rel(c, c_expected, flteps, "srotg(case 208)");
  194.    gsl_test_rel(s, s_expected, flteps, "srotg(case 209)");
  195.   };
  196.  
  197.  
  198.   {
  199.    float a = -1.0f;
  200.    float b = 0.1f;
  201.    float c;
  202.    float s;
  203.    float r_expected = -1.00498756211f;
  204.    float z_expected = -0.099503719021f;
  205.    float c_expected = 0.99503719021f;
  206.    float s_expected = -0.099503719021f;
  207.    cblas_srotg(&a, &b, &c, &s);
  208.    gsl_test_rel(a, r_expected, flteps, "srotg(case 210)");
  209.    gsl_test_rel(b, z_expected, flteps, "srotg(case 211)");
  210.    gsl_test_rel(c, c_expected, flteps, "srotg(case 212)");
  211.    gsl_test_rel(s, s_expected, flteps, "srotg(case 213)");
  212.   };
  213.  
  214.  
  215.   {
  216.    float a = -1.0f;
  217.    float b = 1.0f;
  218.    float c;
  219.    float s;
  220.    float r_expected = 1.41421356237f;
  221.    float z_expected = -1.41421356237f;
  222.    float c_expected = -0.707106781187f;
  223.    float s_expected = 0.707106781187f;
  224.    cblas_srotg(&a, &b, &c, &s);
  225.    gsl_test_rel(a, r_expected, flteps, "srotg(case 214)");
  226.    gsl_test_rel(b, z_expected, flteps, "srotg(case 215)");
  227.    gsl_test_rel(c, c_expected, flteps, "srotg(case 216)");
  228.    gsl_test_rel(s, s_expected, flteps, "srotg(case 217)");
  229.   };
  230.  
  231.  
  232.   {
  233.    float a = -1.0f;
  234.    float b = 1.5f;
  235.    float c;
  236.    float s;
  237.    float r_expected = 1.80277563773f;
  238.    float z_expected = -1.80277563773f;
  239.    float c_expected = -0.554700196225f;
  240.    float s_expected = 0.832050294338f;
  241.    cblas_srotg(&a, &b, &c, &s);
  242.    gsl_test_rel(a, r_expected, flteps, "srotg(case 218)");
  243.    gsl_test_rel(b, z_expected, flteps, "srotg(case 219)");
  244.    gsl_test_rel(c, c_expected, flteps, "srotg(case 220)");
  245.    gsl_test_rel(s, s_expected, flteps, "srotg(case 221)");
  246.   };
  247.  
  248.  
  249.   {
  250.    float a = -0.1f;
  251.    float b = -1.5f;
  252.    float c;
  253.    float s;
  254.    float r_expected = -1.50332963784f;
  255.    float z_expected = 15.0332963784f;
  256.    float c_expected = 0.0665190105238f;
  257.    float s_expected = 0.997785157857f;
  258.    cblas_srotg(&a, &b, &c, &s);
  259.    gsl_test_rel(a, r_expected, flteps, "srotg(case 222)");
  260.    gsl_test_rel(b, z_expected, flteps, "srotg(case 223)");
  261.    gsl_test_rel(c, c_expected, flteps, "srotg(case 224)");
  262.    gsl_test_rel(s, s_expected, flteps, "srotg(case 225)");
  263.   };
  264.  
  265.  
  266.   {
  267.    float a = -0.1f;
  268.    float b = -1.0f;
  269.    float c;
  270.    float s;
  271.    float r_expected = -1.00498756211f;
  272.    float z_expected = 10.0498756211f;
  273.    float c_expected = 0.099503719021f;
  274.    float s_expected = 0.99503719021f;
  275.    cblas_srotg(&a, &b, &c, &s);
  276.    gsl_test_rel(a, r_expected, flteps, "srotg(case 226)");
  277.    gsl_test_rel(b, z_expected, flteps, "srotg(case 227)");
  278.    gsl_test_rel(c, c_expected, flteps, "srotg(case 228)");
  279.    gsl_test_rel(s, s_expected, flteps, "srotg(case 229)");
  280.   };
  281.  
  282.  
  283.   {
  284.    float a = -0.1f;
  285.    float b = -0.1f;
  286.    float c;
  287.    float s;
  288.    float r_expected = -0.141421356237f;
  289.    float z_expected = 1.41421356237f;
  290.    float c_expected = 0.707106781187f;
  291.    float s_expected = 0.707106781187f;
  292.    cblas_srotg(&a, &b, &c, &s);
  293.    gsl_test_rel(a, r_expected, flteps, "srotg(case 230)");
  294.    gsl_test_rel(b, z_expected, flteps, "srotg(case 231)");
  295.    gsl_test_rel(c, c_expected, flteps, "srotg(case 232)");
  296.    gsl_test_rel(s, s_expected, flteps, "srotg(case 233)");
  297.   };
  298.  
  299.  
  300.   {
  301.    float a = -0.1f;
  302.    float b = 0.0f;
  303.    float c;
  304.    float s;
  305.    float r_expected = -0.1f;
  306.    float z_expected = -0.0f;
  307.    float c_expected = 1.0f;
  308.    float s_expected = -0.0f;
  309.    cblas_srotg(&a, &b, &c, &s);
  310.    gsl_test_rel(a, r_expected, flteps, "srotg(case 234)");
  311.    gsl_test_rel(b, z_expected, flteps, "srotg(case 235)");
  312.    gsl_test_rel(c, c_expected, flteps, "srotg(case 236)");
  313.    gsl_test_rel(s, s_expected, flteps, "srotg(case 237)");
  314.   };
  315.  
  316.  
  317.   {
  318.    float a = -0.1f;
  319.    float b = 0.1f;
  320.    float c;
  321.    float s;
  322.    float r_expected = 0.141421356237f;
  323.    float z_expected = -1.41421356237f;
  324.    float c_expected = -0.707106781187f;
  325.    float s_expected = 0.707106781187f;
  326.    cblas_srotg(&a, &b, &c, &s);
  327.    gsl_test_rel(a, r_expected, flteps, "srotg(case 238)");
  328.    gsl_test_rel(b, z_expected, flteps, "srotg(case 239)");
  329.    gsl_test_rel(c, c_expected, flteps, "srotg(case 240)");
  330.    gsl_test_rel(s, s_expected, flteps, "srotg(case 241)");
  331.   };
  332.  
  333.  
  334.   {
  335.    float a = -0.1f;
  336.    float b = 1.0f;
  337.    float c;
  338.    float s;
  339.    float r_expected = 1.00498756211f;
  340.    float z_expected = -10.0498756211f;
  341.    float c_expected = -0.099503719021f;
  342.    float s_expected = 0.99503719021f;
  343.    cblas_srotg(&a, &b, &c, &s);
  344.    gsl_test_rel(a, r_expected, flteps, "srotg(case 242)");
  345.    gsl_test_rel(b, z_expected, flteps, "srotg(case 243)");
  346.    gsl_test_rel(c, c_expected, flteps, "srotg(case 244)");
  347.    gsl_test_rel(s, s_expected, flteps, "srotg(case 245)");
  348.   };
  349.  
  350.  
  351.   {
  352.    float a = -0.1f;
  353.    float b = 1.5f;
  354.    float c;
  355.    float s;
  356.    float r_expected = 1.50332963784f;
  357.    float z_expected = -15.0332963784f;
  358.    float c_expected = -0.0665190105238f;
  359.    float s_expected = 0.997785157857f;
  360.    cblas_srotg(&a, &b, &c, &s);
  361.    gsl_test_rel(a, r_expected, flteps, "srotg(case 246)");
  362.    gsl_test_rel(b, z_expected, flteps, "srotg(case 247)");
  363.    gsl_test_rel(c, c_expected, flteps, "srotg(case 248)");
  364.    gsl_test_rel(s, s_expected, flteps, "srotg(case 249)");
  365.   };
  366.  
  367.  
  368.   {
  369.    float a = 0.0f;
  370.    float b = -1.5f;
  371.    float c;
  372.    float s;
  373.    float r_expected = -1.5f;
  374.    float z_expected = 1.0f;
  375.    float c_expected = -0.0f;
  376.    float s_expected = 1.0f;
  377.    cblas_srotg(&a, &b, &c, &s);
  378.    gsl_test_rel(a, r_expected, flteps, "srotg(case 250)");
  379.    gsl_test_rel(b, z_expected, flteps, "srotg(case 251)");
  380.    gsl_test_rel(c, c_expected, flteps, "srotg(case 252)");
  381.    gsl_test_rel(s, s_expected, flteps, "srotg(case 253)");
  382.   };
  383.  
  384.  
  385.   {
  386.    float a = 0.0f;
  387.    float b = -1.0f;
  388.    float c;
  389.    float s;
  390.    float r_expected = -1.0f;
  391.    float z_expected = 1.0f;
  392.    float c_expected = -0.0f;
  393.    float s_expected = 1.0f;
  394.    cblas_srotg(&a, &b, &c, &s);
  395.    gsl_test_rel(a, r_expected, flteps, "srotg(case 254)");
  396.    gsl_test_rel(b, z_expected, flteps, "srotg(case 255)");
  397.    gsl_test_rel(c, c_expected, flteps, "srotg(case 256)");
  398.    gsl_test_rel(s, s_expected, flteps, "srotg(case 257)");
  399.   };
  400.  
  401.  
  402.   {
  403.    float a = 0.0f;
  404.    float b = -0.1f;
  405.    float c;
  406.    float s;
  407.    float r_expected = -0.1f;
  408.    float z_expected = 1.0f;
  409.    float c_expected = -0.0f;
  410.    float s_expected = 1.0f;
  411.    cblas_srotg(&a, &b, &c, &s);
  412.    gsl_test_rel(a, r_expected, flteps, "srotg(case 258)");
  413.    gsl_test_rel(b, z_expected, flteps, "srotg(case 259)");
  414.    gsl_test_rel(c, c_expected, flteps, "srotg(case 260)");
  415.    gsl_test_rel(s, s_expected, flteps, "srotg(case 261)");
  416.   };
  417.  
  418.  
  419.   {
  420.    float a = 0.0f;
  421.    float b = 0.0f;
  422.    float c;
  423.    float s;
  424.    float r_expected = 0.0f;
  425.    float z_expected = 0.0f;
  426.    float c_expected = 1.0f;
  427.    float s_expected = 0.0f;
  428.    cblas_srotg(&a, &b, &c, &s);
  429.    gsl_test_rel(a, r_expected, flteps, "srotg(case 262)");
  430.    gsl_test_rel(b, z_expected, flteps, "srotg(case 263)");
  431.    gsl_test_rel(c, c_expected, flteps, "srotg(case 264)");
  432.    gsl_test_rel(s, s_expected, flteps, "srotg(case 265)");
  433.   };
  434.  
  435.  
  436.   {
  437.    float a = 0.0f;
  438.    float b = 0.1f;
  439.    float c;
  440.    float s;
  441.    float r_expected = 0.1f;
  442.    float z_expected = 1.0f;
  443.    float c_expected = 0.0f;
  444.    float s_expected = 1.0f;
  445.    cblas_srotg(&a, &b, &c, &s);
  446.    gsl_test_rel(a, r_expected, flteps, "srotg(case 266)");
  447.    gsl_test_rel(b, z_expected, flteps, "srotg(case 267)");
  448.    gsl_test_rel(c, c_expected, flteps, "srotg(case 268)");
  449.    gsl_test_rel(s, s_expected, flteps, "srotg(case 269)");
  450.   };
  451.  
  452.  
  453.   {
  454.    float a = 0.0f;
  455.    float b = 1.0f;
  456.    float c;
  457.    float s;
  458.    float r_expected = 1.0f;
  459.    float z_expected = 1.0f;
  460.    float c_expected = 0.0f;
  461.    float s_expected = 1.0f;
  462.    cblas_srotg(&a, &b, &c, &s);
  463.    gsl_test_rel(a, r_expected, flteps, "srotg(case 270)");
  464.    gsl_test_rel(b, z_expected, flteps, "srotg(case 271)");
  465.    gsl_test_rel(c, c_expected, flteps, "srotg(case 272)");
  466.    gsl_test_rel(s, s_expected, flteps, "srotg(case 273)");
  467.   };
  468.  
  469.  
  470.   {
  471.    float a = 0.0f;
  472.    float b = 1.5f;
  473.    float c;
  474.    float s;
  475.    float r_expected = 1.5f;
  476.    float z_expected = 1.0f;
  477.    float c_expected = 0.0f;
  478.    float s_expected = 1.0f;
  479.    cblas_srotg(&a, &b, &c, &s);
  480.    gsl_test_rel(a, r_expected, flteps, "srotg(case 274)");
  481.    gsl_test_rel(b, z_expected, flteps, "srotg(case 275)");
  482.    gsl_test_rel(c, c_expected, flteps, "srotg(case 276)");
  483.    gsl_test_rel(s, s_expected, flteps, "srotg(case 277)");
  484.   };
  485.  
  486.  
  487.   {
  488.    float a = 0.1f;
  489.    float b = -1.5f;
  490.    float c;
  491.    float s;
  492.    float r_expected = -1.50332963784f;
  493.    float z_expected = -15.0332963784f;
  494.    float c_expected = -0.0665190105238f;
  495.    float s_expected = 0.997785157857f;
  496.    cblas_srotg(&a, &b, &c, &s);
  497.    gsl_test_rel(a, r_expected, flteps, "srotg(case 278)");
  498.    gsl_test_rel(b, z_expected, flteps, "srotg(case 279)");
  499.    gsl_test_rel(c, c_expected, flteps, "srotg(case 280)");
  500.    gsl_test_rel(s, s_expected, flteps, "srotg(case 281)");
  501.   };
  502.  
  503.  
  504.   {
  505.    float a = 0.1f;
  506.    float b = -1.0f;
  507.    float c;
  508.    float s;
  509.    float r_expected = -1.00498756211f;
  510.    float z_expected = -10.0498756211f;
  511.    float c_expected = -0.099503719021f;
  512.    float s_expected = 0.99503719021f;
  513.    cblas_srotg(&a, &b, &c, &s);
  514.    gsl_test_rel(a, r_expected, flteps, "srotg(case 282)");
  515.    gsl_test_rel(b, z_expected, flteps, "srotg(case 283)");
  516.    gsl_test_rel(c, c_expected, flteps, "srotg(case 284)");
  517.    gsl_test_rel(s, s_expected, flteps, "srotg(case 285)");
  518.   };
  519.  
  520.  
  521.   {
  522.    float a = 0.1f;
  523.    float b = -0.1f;
  524.    float c;
  525.    float s;
  526.    float r_expected = -0.141421356237f;
  527.    float z_expected = -1.41421356237f;
  528.    float c_expected = -0.707106781187f;
  529.    float s_expected = 0.707106781187f;
  530.    cblas_srotg(&a, &b, &c, &s);
  531.    gsl_test_rel(a, r_expected, flteps, "srotg(case 286)");
  532.    gsl_test_rel(b, z_expected, flteps, "srotg(case 287)");
  533.    gsl_test_rel(c, c_expected, flteps, "srotg(case 288)");
  534.    gsl_test_rel(s, s_expected, flteps, "srotg(case 289)");
  535.   };
  536.  
  537.  
  538.   {
  539.    float a = 0.1f;
  540.    float b = 0.0f;
  541.    float c;
  542.    float s;
  543.    float r_expected = 0.1f;
  544.    float z_expected = 0.0f;
  545.    float c_expected = 1.0f;
  546.    float s_expected = 0.0f;
  547.    cblas_srotg(&a, &b, &c, &s);
  548.    gsl_test_rel(a, r_expected, flteps, "srotg(case 290)");
  549.    gsl_test_rel(b, z_expected, flteps, "srotg(case 291)");
  550.    gsl_test_rel(c, c_expected, flteps, "srotg(case 292)");
  551.    gsl_test_rel(s, s_expected, flteps, "srotg(case 293)");
  552.   };
  553.  
  554.  
  555.   {
  556.    float a = 0.1f;
  557.    float b = 0.1f;
  558.    float c;
  559.    float s;
  560.    float r_expected = 0.141421356237f;
  561.    float z_expected = 1.41421356237f;
  562.    float c_expected = 0.707106781187f;
  563.    float s_expected = 0.707106781187f;
  564.    cblas_srotg(&a, &b, &c, &s);
  565.    gsl_test_rel(a, r_expected, flteps, "srotg(case 294)");
  566.    gsl_test_rel(b, z_expected, flteps, "srotg(case 295)");
  567.    gsl_test_rel(c, c_expected, flteps, "srotg(case 296)");
  568.    gsl_test_rel(s, s_expected, flteps, "srotg(case 297)");
  569.   };
  570.  
  571.  
  572.   {
  573.    float a = 0.1f;
  574.    float b = 1.0f;
  575.    float c;
  576.    float s;
  577.    float r_expected = 1.00498756211f;
  578.    float z_expected = 10.0498756211f;
  579.    float c_expected = 0.099503719021f;
  580.    float s_expected = 0.99503719021f;
  581.    cblas_srotg(&a, &b, &c, &s);
  582.    gsl_test_rel(a, r_expected, flteps, "srotg(case 298)");
  583.    gsl_test_rel(b, z_expected, flteps, "srotg(case 299)");
  584.    gsl_test_rel(c, c_expected, flteps, "srotg(case 300)");
  585.    gsl_test_rel(s, s_expected, flteps, "srotg(case 301)");
  586.   };
  587.  
  588.  
  589.   {
  590.    float a = 0.1f;
  591.    float b = 1.5f;
  592.    float c;
  593.    float s;
  594.    float r_expected = 1.50332963784f;
  595.    float z_expected = 15.0332963784f;
  596.    float c_expected = 0.0665190105238f;
  597.    float s_expected = 0.997785157857f;
  598.    cblas_srotg(&a, &b, &c, &s);
  599.    gsl_test_rel(a, r_expected, flteps, "srotg(case 302)");
  600.    gsl_test_rel(b, z_expected, flteps, "srotg(case 303)");
  601.    gsl_test_rel(c, c_expected, flteps, "srotg(case 304)");
  602.    gsl_test_rel(s, s_expected, flteps, "srotg(case 305)");
  603.   };
  604.  
  605.  
  606.   {
  607.    float a = 1.0f;
  608.    float b = -1.5f;
  609.    float c;
  610.    float s;
  611.    float r_expected = -1.80277563773f;
  612.    float z_expected = -1.80277563773f;
  613.    float c_expected = -0.554700196225f;
  614.    float s_expected = 0.832050294338f;
  615.    cblas_srotg(&a, &b, &c, &s);
  616.    gsl_test_rel(a, r_expected, flteps, "srotg(case 306)");
  617.    gsl_test_rel(b, z_expected, flteps, "srotg(case 307)");
  618.    gsl_test_rel(c, c_expected, flteps, "srotg(case 308)");
  619.    gsl_test_rel(s, s_expected, flteps, "srotg(case 309)");
  620.   };
  621.  
  622.  
  623.   {
  624.    float a = 1.0f;
  625.    float b = -1.0f;
  626.    float c;
  627.    float s;
  628.    float r_expected = -1.41421356237f;
  629.    float z_expected = -1.41421356237f;
  630.    float c_expected = -0.707106781187f;
  631.    float s_expected = 0.707106781187f;
  632.    cblas_srotg(&a, &b, &c, &s);
  633.    gsl_test_rel(a, r_expected, flteps, "srotg(case 310)");
  634.    gsl_test_rel(b, z_expected, flteps, "srotg(case 311)");
  635.    gsl_test_rel(c, c_expected, flteps, "srotg(case 312)");
  636.    gsl_test_rel(s, s_expected, flteps, "srotg(case 313)");
  637.   };
  638.  
  639.  
  640.   {
  641.    float a = 1.0f;
  642.    float b = -0.1f;
  643.    float c;
  644.    float s;
  645.    float r_expected = 1.00498756211f;
  646.    float z_expected = -0.099503719021f;
  647.    float c_expected = 0.99503719021f;
  648.    float s_expected = -0.099503719021f;
  649.    cblas_srotg(&a, &b, &c, &s);
  650.    gsl_test_rel(a, r_expected, flteps, "srotg(case 314)");
  651.    gsl_test_rel(b, z_expected, flteps, "srotg(case 315)");
  652.    gsl_test_rel(c, c_expected, flteps, "srotg(case 316)");
  653.    gsl_test_rel(s, s_expected, flteps, "srotg(case 317)");
  654.   };
  655.  
  656.  
  657.   {
  658.    float a = 1.0f;
  659.    float b = 0.0f;
  660.    float c;
  661.    float s;
  662.    float r_expected = 1.0f;
  663.    float z_expected = 0.0f;
  664.    float c_expected = 1.0f;
  665.    float s_expected = 0.0f;
  666.    cblas_srotg(&a, &b, &c, &s);
  667.    gsl_test_rel(a, r_expected, flteps, "srotg(case 318)");
  668.    gsl_test_rel(b, z_expected, flteps, "srotg(case 319)");
  669.    gsl_test_rel(c, c_expected, flteps, "srotg(case 320)");
  670.    gsl_test_rel(s, s_expected, flteps, "srotg(case 321)");
  671.   };
  672.  
  673.  
  674.   {
  675.    float a = 1.0f;
  676.    float b = 0.1f;
  677.    float c;
  678.    float s;
  679.    float r_expected = 1.00498756211f;
  680.    float z_expected = 0.099503719021f;
  681.    float c_expected = 0.99503719021f;
  682.    float s_expected = 0.099503719021f;
  683.    cblas_srotg(&a, &b, &c, &s);
  684.    gsl_test_rel(a, r_expected, flteps, "srotg(case 322)");
  685.    gsl_test_rel(b, z_expected, flteps, "srotg(case 323)");
  686.    gsl_test_rel(c, c_expected, flteps, "srotg(case 324)");
  687.    gsl_test_rel(s, s_expected, flteps, "srotg(case 325)");
  688.   };
  689.  
  690.  
  691.   {
  692.    float a = 1.0f;
  693.    float b = 1.0f;
  694.    float c;
  695.    float s;
  696.    float r_expected = 1.41421356237f;
  697.    float z_expected = 1.41421356237f;
  698.    float c_expected = 0.707106781187f;
  699.    float s_expected = 0.707106781187f;
  700.    cblas_srotg(&a, &b, &c, &s);
  701.    gsl_test_rel(a, r_expected, flteps, "srotg(case 326)");
  702.    gsl_test_rel(b, z_expected, flteps, "srotg(case 327)");
  703.    gsl_test_rel(c, c_expected, flteps, "srotg(case 328)");
  704.    gsl_test_rel(s, s_expected, flteps, "srotg(case 329)");
  705.   };
  706.  
  707.  
  708.   {
  709.    float a = 1.0f;
  710.    float b = 1.5f;
  711.    float c;
  712.    float s;
  713.    float r_expected = 1.80277563773f;
  714.    float z_expected = 1.80277563773f;
  715.    float c_expected = 0.554700196225f;
  716.    float s_expected = 0.832050294338f;
  717.    cblas_srotg(&a, &b, &c, &s);
  718.    gsl_test_rel(a, r_expected, flteps, "srotg(case 330)");
  719.    gsl_test_rel(b, z_expected, flteps, "srotg(case 331)");
  720.    gsl_test_rel(c, c_expected, flteps, "srotg(case 332)");
  721.    gsl_test_rel(s, s_expected, flteps, "srotg(case 333)");
  722.   };
  723.  
  724.  
  725.   {
  726.    float a = 1.5f;
  727.    float b = -1.5f;
  728.    float c;
  729.    float s;
  730.    float r_expected = -2.12132034356f;
  731.    float z_expected = -1.41421356237f;
  732.    float c_expected = -0.707106781187f;
  733.    float s_expected = 0.707106781187f;
  734.    cblas_srotg(&a, &b, &c, &s);
  735.    gsl_test_rel(a, r_expected, flteps, "srotg(case 334)");
  736.    gsl_test_rel(b, z_expected, flteps, "srotg(case 335)");
  737.    gsl_test_rel(c, c_expected, flteps, "srotg(case 336)");
  738.    gsl_test_rel(s, s_expected, flteps, "srotg(case 337)");
  739.   };
  740.  
  741.  
  742.   {
  743.    float a = 1.5f;
  744.    float b = -1.0f;
  745.    float c;
  746.    float s;
  747.    float r_expected = 1.80277563773f;
  748.    float z_expected = -0.554700196225f;
  749.    float c_expected = 0.832050294338f;
  750.    float s_expected = -0.554700196225f;
  751.    cblas_srotg(&a, &b, &c, &s);
  752.    gsl_test_rel(a, r_expected, flteps, "srotg(case 338)");
  753.    gsl_test_rel(b, z_expected, flteps, "srotg(case 339)");
  754.    gsl_test_rel(c, c_expected, flteps, "srotg(case 340)");
  755.    gsl_test_rel(s, s_expected, flteps, "srotg(case 341)");
  756.   };
  757.  
  758.  
  759.   {
  760.    float a = 1.5f;
  761.    float b = -0.1f;
  762.    float c;
  763.    float s;
  764.    float r_expected = 1.50332963784f;
  765.    float z_expected = -0.0665190105238f;
  766.    float c_expected = 0.997785157857f;
  767.    float s_expected = -0.0665190105238f;
  768.    cblas_srotg(&a, &b, &c, &s);
  769.    gsl_test_rel(a, r_expected, flteps, "srotg(case 342)");
  770.    gsl_test_rel(b, z_expected, flteps, "srotg(case 343)");
  771.    gsl_test_rel(c, c_expected, flteps, "srotg(case 344)");
  772.    gsl_test_rel(s, s_expected, flteps, "srotg(case 345)");
  773.   };
  774.  
  775.  
  776.   {
  777.    float a = 1.5f;
  778.    float b = 0.0f;
  779.    float c;
  780.    float s;
  781.    float r_expected = 1.5f;
  782.    float z_expected = 0.0f;
  783.    float c_expected = 1.0f;
  784.    float s_expected = 0.0f;
  785.    cblas_srotg(&a, &b, &c, &s);
  786.    gsl_test_rel(a, r_expected, flteps, "srotg(case 346)");
  787.    gsl_test_rel(b, z_expected, flteps, "srotg(case 347)");
  788.    gsl_test_rel(c, c_expected, flteps, "srotg(case 348)");
  789.    gsl_test_rel(s, s_expected, flteps, "srotg(case 349)");
  790.   };
  791.  
  792.  
  793.   {
  794.    float a = 1.5f;
  795.    float b = 0.1f;
  796.    float c;
  797.    float s;
  798.    float r_expected = 1.50332963784f;
  799.    float z_expected = 0.0665190105238f;
  800.    float c_expected = 0.997785157857f;
  801.    float s_expected = 0.0665190105238f;
  802.    cblas_srotg(&a, &b, &c, &s);
  803.    gsl_test_rel(a, r_expected, flteps, "srotg(case 350)");
  804.    gsl_test_rel(b, z_expected, flteps, "srotg(case 351)");
  805.    gsl_test_rel(c, c_expected, flteps, "srotg(case 352)");
  806.    gsl_test_rel(s, s_expected, flteps, "srotg(case 353)");
  807.   };
  808.  
  809.  
  810.   {
  811.    float a = 1.5f;
  812.    float b = 1.0f;
  813.    float c;
  814.    float s;
  815.    float r_expected = 1.80277563773f;
  816.    float z_expected = 0.554700196225f;
  817.    float c_expected = 0.832050294338f;
  818.    float s_expected = 0.554700196225f;
  819.    cblas_srotg(&a, &b, &c, &s);
  820.    gsl_test_rel(a, r_expected, flteps, "srotg(case 354)");
  821.    gsl_test_rel(b, z_expected, flteps, "srotg(case 355)");
  822.    gsl_test_rel(c, c_expected, flteps, "srotg(case 356)");
  823.    gsl_test_rel(s, s_expected, flteps, "srotg(case 357)");
  824.   };
  825.  
  826.  
  827.   {
  828.    float a = 1.5f;
  829.    float b = 1.5f;
  830.    float c;
  831.    float s;
  832.    float r_expected = 2.12132034356f;
  833.    float z_expected = 1.41421356237f;
  834.    float c_expected = 0.707106781187f;
  835.    float s_expected = 0.707106781187f;
  836.    cblas_srotg(&a, &b, &c, &s);
  837.    gsl_test_rel(a, r_expected, flteps, "srotg(case 358)");
  838.    gsl_test_rel(b, z_expected, flteps, "srotg(case 359)");
  839.    gsl_test_rel(c, c_expected, flteps, "srotg(case 360)");
  840.    gsl_test_rel(s, s_expected, flteps, "srotg(case 361)");
  841.   };
  842.  
  843.  
  844.   {
  845.    double a = -1.5;
  846.    double b = -1.5;
  847.    double c;
  848.    double s;
  849.    double r_expected = -2.12132034356;
  850.    double z_expected = 1.41421356237;
  851.    double c_expected = 0.707106781187;
  852.    double s_expected = 0.707106781187;
  853.    cblas_drotg(&a, &b, &c, &s);
  854.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 362)");
  855.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 363)");
  856.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 364)");
  857.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 365)");
  858.   };
  859.  
  860.  
  861.   {
  862.    double a = -1.5;
  863.    double b = -1;
  864.    double c;
  865.    double s;
  866.    double r_expected = -1.80277563773;
  867.    double z_expected = 0.554700196225;
  868.    double c_expected = 0.832050294338;
  869.    double s_expected = 0.554700196225;
  870.    cblas_drotg(&a, &b, &c, &s);
  871.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 366)");
  872.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 367)");
  873.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 368)");
  874.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 369)");
  875.   };
  876.  
  877.  
  878.   {
  879.    double a = -1.5;
  880.    double b = -0.1;
  881.    double c;
  882.    double s;
  883.    double r_expected = -1.50332963784;
  884.    double z_expected = 0.0665190105238;
  885.    double c_expected = 0.997785157857;
  886.    double s_expected = 0.0665190105238;
  887.    cblas_drotg(&a, &b, &c, &s);
  888.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 370)");
  889.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 371)");
  890.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 372)");
  891.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 373)");
  892.   };
  893.  
  894.  
  895.   {
  896.    double a = -1.5;
  897.    double b = 0;
  898.    double c;
  899.    double s;
  900.    double r_expected = -1.5;
  901.    double z_expected = -0;
  902.    double c_expected = 1;
  903.    double s_expected = -0;
  904.    cblas_drotg(&a, &b, &c, &s);
  905.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 374)");
  906.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 375)");
  907.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 376)");
  908.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 377)");
  909.   };
  910.  
  911.  
  912.   {
  913.    double a = -1.5;
  914.    double b = 0.1;
  915.    double c;
  916.    double s;
  917.    double r_expected = -1.50332963784;
  918.    double z_expected = -0.0665190105238;
  919.    double c_expected = 0.997785157857;
  920.    double s_expected = -0.0665190105238;
  921.    cblas_drotg(&a, &b, &c, &s);
  922.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 378)");
  923.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 379)");
  924.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 380)");
  925.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 381)");
  926.   };
  927.  
  928.  
  929.   {
  930.    double a = -1.5;
  931.    double b = 1;
  932.    double c;
  933.    double s;
  934.    double r_expected = -1.80277563773;
  935.    double z_expected = -0.554700196225;
  936.    double c_expected = 0.832050294338;
  937.    double s_expected = -0.554700196225;
  938.    cblas_drotg(&a, &b, &c, &s);
  939.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 382)");
  940.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 383)");
  941.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 384)");
  942.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 385)");
  943.   };
  944.  
  945.  
  946.   {
  947.    double a = -1.5;
  948.    double b = 1.5;
  949.    double c;
  950.    double s;
  951.    double r_expected = 2.12132034356;
  952.    double z_expected = -1.41421356237;
  953.    double c_expected = -0.707106781187;
  954.    double s_expected = 0.707106781187;
  955.    cblas_drotg(&a, &b, &c, &s);
  956.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 386)");
  957.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 387)");
  958.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 388)");
  959.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 389)");
  960.   };
  961.  
  962.  
  963.   {
  964.    double a = -1;
  965.    double b = -1.5;
  966.    double c;
  967.    double s;
  968.    double r_expected = -1.80277563773;
  969.    double z_expected = 1.80277563773;
  970.    double c_expected = 0.554700196225;
  971.    double s_expected = 0.832050294338;
  972.    cblas_drotg(&a, &b, &c, &s);
  973.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 390)");
  974.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 391)");
  975.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 392)");
  976.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 393)");
  977.   };
  978.  
  979.  
  980.   {
  981.    double a = -1;
  982.    double b = -1;
  983.    double c;
  984.    double s;
  985.    double r_expected = -1.41421356237;
  986.    double z_expected = 1.41421356237;
  987.    double c_expected = 0.707106781187;
  988.    double s_expected = 0.707106781187;
  989.    cblas_drotg(&a, &b, &c, &s);
  990.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 394)");
  991.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 395)");
  992.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 396)");
  993.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 397)");
  994.   };
  995.  
  996.  
  997.   {
  998.    double a = -1;
  999.    double b = -0.1;
  1000.    double c;
  1001.    double s;
  1002.    double r_expected = -1.00498756211;
  1003.    double z_expected = 0.099503719021;
  1004.    double c_expected = 0.99503719021;
  1005.    double s_expected = 0.099503719021;
  1006.    cblas_drotg(&a, &b, &c, &s);
  1007.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 398)");
  1008.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 399)");
  1009.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 400)");
  1010.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 401)");
  1011.   };
  1012.  
  1013.  
  1014.   {
  1015.    double a = -1;
  1016.    double b = 0;
  1017.    double c;
  1018.    double s;
  1019.    double r_expected = -1;
  1020.    double z_expected = -0;
  1021.    double c_expected = 1;
  1022.    double s_expected = -0;
  1023.    cblas_drotg(&a, &b, &c, &s);
  1024.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 402)");
  1025.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 403)");
  1026.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 404)");
  1027.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 405)");
  1028.   };
  1029.  
  1030.  
  1031.   {
  1032.    double a = -1;
  1033.    double b = 0.1;
  1034.    double c;
  1035.    double s;
  1036.    double r_expected = -1.00498756211;
  1037.    double z_expected = -0.099503719021;
  1038.    double c_expected = 0.99503719021;
  1039.    double s_expected = -0.099503719021;
  1040.    cblas_drotg(&a, &b, &c, &s);
  1041.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 406)");
  1042.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 407)");
  1043.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 408)");
  1044.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 409)");
  1045.   };
  1046.  
  1047.  
  1048.   {
  1049.    double a = -1;
  1050.    double b = 1;
  1051.    double c;
  1052.    double s;
  1053.    double r_expected = 1.41421356237;
  1054.    double z_expected = -1.41421356237;
  1055.    double c_expected = -0.707106781187;
  1056.    double s_expected = 0.707106781187;
  1057.    cblas_drotg(&a, &b, &c, &s);
  1058.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 410)");
  1059.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 411)");
  1060.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 412)");
  1061.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 413)");
  1062.   };
  1063.  
  1064.  
  1065.   {
  1066.    double a = -1;
  1067.    double b = 1.5;
  1068.    double c;
  1069.    double s;
  1070.    double r_expected = 1.80277563773;
  1071.    double z_expected = -1.80277563773;
  1072.    double c_expected = -0.554700196225;
  1073.    double s_expected = 0.832050294338;
  1074.    cblas_drotg(&a, &b, &c, &s);
  1075.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 414)");
  1076.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 415)");
  1077.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 416)");
  1078.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 417)");
  1079.   };
  1080.  
  1081.  
  1082.   {
  1083.    double a = -0.1;
  1084.    double b = -1.5;
  1085.    double c;
  1086.    double s;
  1087.    double r_expected = -1.50332963784;
  1088.    double z_expected = 15.0332963784;
  1089.    double c_expected = 0.0665190105238;
  1090.    double s_expected = 0.997785157857;
  1091.    cblas_drotg(&a, &b, &c, &s);
  1092.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 418)");
  1093.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 419)");
  1094.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 420)");
  1095.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 421)");
  1096.   };
  1097.  
  1098.  
  1099.   {
  1100.    double a = -0.1;
  1101.    double b = -1;
  1102.    double c;
  1103.    double s;
  1104.    double r_expected = -1.00498756211;
  1105.    double z_expected = 10.0498756211;
  1106.    double c_expected = 0.099503719021;
  1107.    double s_expected = 0.99503719021;
  1108.    cblas_drotg(&a, &b, &c, &s);
  1109.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 422)");
  1110.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 423)");
  1111.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 424)");
  1112.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 425)");
  1113.   };
  1114.  
  1115.  
  1116.   {
  1117.    double a = -0.1;
  1118.    double b = -0.1;
  1119.    double c;
  1120.    double s;
  1121.    double r_expected = -0.141421356237;
  1122.    double z_expected = 1.41421356237;
  1123.    double c_expected = 0.707106781187;
  1124.    double s_expected = 0.707106781187;
  1125.    cblas_drotg(&a, &b, &c, &s);
  1126.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 426)");
  1127.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 427)");
  1128.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 428)");
  1129.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 429)");
  1130.   };
  1131.  
  1132.  
  1133.   {
  1134.    double a = -0.1;
  1135.    double b = 0;
  1136.    double c;
  1137.    double s;
  1138.    double r_expected = -0.1;
  1139.    double z_expected = -0;
  1140.    double c_expected = 1;
  1141.    double s_expected = -0;
  1142.    cblas_drotg(&a, &b, &c, &s);
  1143.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 430)");
  1144.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 431)");
  1145.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 432)");
  1146.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 433)");
  1147.   };
  1148.  
  1149.  
  1150.   {
  1151.    double a = -0.1;
  1152.    double b = 0.1;
  1153.    double c;
  1154.    double s;
  1155.    double r_expected = 0.141421356237;
  1156.    double z_expected = -1.41421356237;
  1157.    double c_expected = -0.707106781187;
  1158.    double s_expected = 0.707106781187;
  1159.    cblas_drotg(&a, &b, &c, &s);
  1160.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 434)");
  1161.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 435)");
  1162.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 436)");
  1163.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 437)");
  1164.   };
  1165.  
  1166.  
  1167.   {
  1168.    double a = -0.1;
  1169.    double b = 1;
  1170.    double c;
  1171.    double s;
  1172.    double r_expected = 1.00498756211;
  1173.    double z_expected = -10.0498756211;
  1174.    double c_expected = -0.099503719021;
  1175.    double s_expected = 0.99503719021;
  1176.    cblas_drotg(&a, &b, &c, &s);
  1177.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 438)");
  1178.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 439)");
  1179.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 440)");
  1180.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 441)");
  1181.   };
  1182.  
  1183.  
  1184.   {
  1185.    double a = -0.1;
  1186.    double b = 1.5;
  1187.    double c;
  1188.    double s;
  1189.    double r_expected = 1.50332963784;
  1190.    double z_expected = -15.0332963784;
  1191.    double c_expected = -0.0665190105238;
  1192.    double s_expected = 0.997785157857;
  1193.    cblas_drotg(&a, &b, &c, &s);
  1194.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 442)");
  1195.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 443)");
  1196.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 444)");
  1197.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 445)");
  1198.   };
  1199.  
  1200.  
  1201.   {
  1202.    double a = 0;
  1203.    double b = -1.5;
  1204.    double c;
  1205.    double s;
  1206.    double r_expected = -1.5;
  1207.    double z_expected = 1;
  1208.    double c_expected = -0;
  1209.    double s_expected = 1;
  1210.    cblas_drotg(&a, &b, &c, &s);
  1211.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 446)");
  1212.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 447)");
  1213.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 448)");
  1214.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 449)");
  1215.   };
  1216.  
  1217.  
  1218.   {
  1219.    double a = 0;
  1220.    double b = -1;
  1221.    double c;
  1222.    double s;
  1223.    double r_expected = -1;
  1224.    double z_expected = 1;
  1225.    double c_expected = -0;
  1226.    double s_expected = 1;
  1227.    cblas_drotg(&a, &b, &c, &s);
  1228.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 450)");
  1229.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 451)");
  1230.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 452)");
  1231.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 453)");
  1232.   };
  1233.  
  1234.  
  1235.   {
  1236.    double a = 0;
  1237.    double b = -0.1;
  1238.    double c;
  1239.    double s;
  1240.    double r_expected = -0.1;
  1241.    double z_expected = 1;
  1242.    double c_expected = -0;
  1243.    double s_expected = 1;
  1244.    cblas_drotg(&a, &b, &c, &s);
  1245.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 454)");
  1246.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 455)");
  1247.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 456)");
  1248.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 457)");
  1249.   };
  1250.  
  1251.  
  1252.   {
  1253.    double a = 0;
  1254.    double b = 0;
  1255.    double c;
  1256.    double s;
  1257.    double r_expected = 0;
  1258.    double z_expected = 0;
  1259.    double c_expected = 1;
  1260.    double s_expected = 0;
  1261.    cblas_drotg(&a, &b, &c, &s);
  1262.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 458)");
  1263.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 459)");
  1264.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 460)");
  1265.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 461)");
  1266.   };
  1267.  
  1268.  
  1269.   {
  1270.    double a = 0;
  1271.    double b = 0.1;
  1272.    double c;
  1273.    double s;
  1274.    double r_expected = 0.1;
  1275.    double z_expected = 1;
  1276.    double c_expected = 0;
  1277.    double s_expected = 1;
  1278.    cblas_drotg(&a, &b, &c, &s);
  1279.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 462)");
  1280.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 463)");
  1281.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 464)");
  1282.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 465)");
  1283.   };
  1284.  
  1285.  
  1286.   {
  1287.    double a = 0;
  1288.    double b = 1;
  1289.    double c;
  1290.    double s;
  1291.    double r_expected = 1;
  1292.    double z_expected = 1;
  1293.    double c_expected = 0;
  1294.    double s_expected = 1;
  1295.    cblas_drotg(&a, &b, &c, &s);
  1296.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 466)");
  1297.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 467)");
  1298.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 468)");
  1299.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 469)");
  1300.   };
  1301.  
  1302.  
  1303.   {
  1304.    double a = 0;
  1305.    double b = 1.5;
  1306.    double c;
  1307.    double s;
  1308.    double r_expected = 1.5;
  1309.    double z_expected = 1;
  1310.    double c_expected = 0;
  1311.    double s_expected = 1;
  1312.    cblas_drotg(&a, &b, &c, &s);
  1313.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 470)");
  1314.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 471)");
  1315.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 472)");
  1316.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 473)");
  1317.   };
  1318.  
  1319.  
  1320.   {
  1321.    double a = 0.1;
  1322.    double b = -1.5;
  1323.    double c;
  1324.    double s;
  1325.    double r_expected = -1.50332963784;
  1326.    double z_expected = -15.0332963784;
  1327.    double c_expected = -0.0665190105238;
  1328.    double s_expected = 0.997785157857;
  1329.    cblas_drotg(&a, &b, &c, &s);
  1330.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 474)");
  1331.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 475)");
  1332.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 476)");
  1333.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 477)");
  1334.   };
  1335.  
  1336.  
  1337.   {
  1338.    double a = 0.1;
  1339.    double b = -1;
  1340.    double c;
  1341.    double s;
  1342.    double r_expected = -1.00498756211;
  1343.    double z_expected = -10.0498756211;
  1344.    double c_expected = -0.099503719021;
  1345.    double s_expected = 0.99503719021;
  1346.    cblas_drotg(&a, &b, &c, &s);
  1347.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 478)");
  1348.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 479)");
  1349.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 480)");
  1350.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 481)");
  1351.   };
  1352.  
  1353.  
  1354.   {
  1355.    double a = 0.1;
  1356.    double b = -0.1;
  1357.    double c;
  1358.    double s;
  1359.    double r_expected = -0.141421356237;
  1360.    double z_expected = -1.41421356237;
  1361.    double c_expected = -0.707106781187;
  1362.    double s_expected = 0.707106781187;
  1363.    cblas_drotg(&a, &b, &c, &s);
  1364.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 482)");
  1365.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 483)");
  1366.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 484)");
  1367.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 485)");
  1368.   };
  1369.  
  1370.  
  1371.   {
  1372.    double a = 0.1;
  1373.    double b = 0;
  1374.    double c;
  1375.    double s;
  1376.    double r_expected = 0.1;
  1377.    double z_expected = 0;
  1378.    double c_expected = 1;
  1379.    double s_expected = 0;
  1380.    cblas_drotg(&a, &b, &c, &s);
  1381.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 486)");
  1382.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 487)");
  1383.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 488)");
  1384.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 489)");
  1385.   };
  1386.  
  1387.  
  1388.   {
  1389.    double a = 0.1;
  1390.    double b = 0.1;
  1391.    double c;
  1392.    double s;
  1393.    double r_expected = 0.141421356237;
  1394.    double z_expected = 1.41421356237;
  1395.    double c_expected = 0.707106781187;
  1396.    double s_expected = 0.707106781187;
  1397.    cblas_drotg(&a, &b, &c, &s);
  1398.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 490)");
  1399.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 491)");
  1400.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 492)");
  1401.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 493)");
  1402.   };
  1403.  
  1404.  
  1405.   {
  1406.    double a = 0.1;
  1407.    double b = 1;
  1408.    double c;
  1409.    double s;
  1410.    double r_expected = 1.00498756211;
  1411.    double z_expected = 10.0498756211;
  1412.    double c_expected = 0.099503719021;
  1413.    double s_expected = 0.99503719021;
  1414.    cblas_drotg(&a, &b, &c, &s);
  1415.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 494)");
  1416.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 495)");
  1417.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 496)");
  1418.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 497)");
  1419.   };
  1420.  
  1421.  
  1422.   {
  1423.    double a = 0.1;
  1424.    double b = 1.5;
  1425.    double c;
  1426.    double s;
  1427.    double r_expected = 1.50332963784;
  1428.    double z_expected = 15.0332963784;
  1429.    double c_expected = 0.0665190105238;
  1430.    double s_expected = 0.997785157857;
  1431.    cblas_drotg(&a, &b, &c, &s);
  1432.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 498)");
  1433.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 499)");
  1434.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 500)");
  1435.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 501)");
  1436.   };
  1437.  
  1438.  
  1439.   {
  1440.    double a = 1;
  1441.    double b = -1.5;
  1442.    double c;
  1443.    double s;
  1444.    double r_expected = -1.80277563773;
  1445.    double z_expected = -1.80277563773;
  1446.    double c_expected = -0.554700196225;
  1447.    double s_expected = 0.832050294338;
  1448.    cblas_drotg(&a, &b, &c, &s);
  1449.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 502)");
  1450.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 503)");
  1451.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 504)");
  1452.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 505)");
  1453.   };
  1454.  
  1455.  
  1456.   {
  1457.    double a = 1;
  1458.    double b = -1;
  1459.    double c;
  1460.    double s;
  1461.    double r_expected = -1.41421356237;
  1462.    double z_expected = -1.41421356237;
  1463.    double c_expected = -0.707106781187;
  1464.    double s_expected = 0.707106781187;
  1465.    cblas_drotg(&a, &b, &c, &s);
  1466.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 506)");
  1467.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 507)");
  1468.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 508)");
  1469.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 509)");
  1470.   };
  1471.  
  1472.  
  1473.   {
  1474.    double a = 1;
  1475.    double b = -0.1;
  1476.    double c;
  1477.    double s;
  1478.    double r_expected = 1.00498756211;
  1479.    double z_expected = -0.099503719021;
  1480.    double c_expected = 0.99503719021;
  1481.    double s_expected = -0.099503719021;
  1482.    cblas_drotg(&a, &b, &c, &s);
  1483.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 510)");
  1484.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 511)");
  1485.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 512)");
  1486.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 513)");
  1487.   };
  1488.  
  1489.  
  1490.   {
  1491.    double a = 1;
  1492.    double b = 0;
  1493.    double c;
  1494.    double s;
  1495.    double r_expected = 1;
  1496.    double z_expected = 0;
  1497.    double c_expected = 1;
  1498.    double s_expected = 0;
  1499.    cblas_drotg(&a, &b, &c, &s);
  1500.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 514)");
  1501.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 515)");
  1502.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 516)");
  1503.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 517)");
  1504.   };
  1505.  
  1506.  
  1507.   {
  1508.    double a = 1;
  1509.    double b = 0.1;
  1510.    double c;
  1511.    double s;
  1512.    double r_expected = 1.00498756211;
  1513.    double z_expected = 0.099503719021;
  1514.    double c_expected = 0.99503719021;
  1515.    double s_expected = 0.099503719021;
  1516.    cblas_drotg(&a, &b, &c, &s);
  1517.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 518)");
  1518.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 519)");
  1519.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 520)");
  1520.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 521)");
  1521.   };
  1522.  
  1523.  
  1524.   {
  1525.    double a = 1;
  1526.    double b = 1;
  1527.    double c;
  1528.    double s;
  1529.    double r_expected = 1.41421356237;
  1530.    double z_expected = 1.41421356237;
  1531.    double c_expected = 0.707106781187;
  1532.    double s_expected = 0.707106781187;
  1533.    cblas_drotg(&a, &b, &c, &s);
  1534.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 522)");
  1535.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 523)");
  1536.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 524)");
  1537.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 525)");
  1538.   };
  1539.  
  1540.  
  1541.   {
  1542.    double a = 1;
  1543.    double b = 1.5;
  1544.    double c;
  1545.    double s;
  1546.    double r_expected = 1.80277563773;
  1547.    double z_expected = 1.80277563773;
  1548.    double c_expected = 0.554700196225;
  1549.    double s_expected = 0.832050294338;
  1550.    cblas_drotg(&a, &b, &c, &s);
  1551.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 526)");
  1552.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 527)");
  1553.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 528)");
  1554.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 529)");
  1555.   };
  1556.  
  1557.  
  1558.   {
  1559.    double a = 1.5;
  1560.    double b = -1.5;
  1561.    double c;
  1562.    double s;
  1563.    double r_expected = -2.12132034356;
  1564.    double z_expected = -1.41421356237;
  1565.    double c_expected = -0.707106781187;
  1566.    double s_expected = 0.707106781187;
  1567.    cblas_drotg(&a, &b, &c, &s);
  1568.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 530)");
  1569.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 531)");
  1570.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 532)");
  1571.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 533)");
  1572.   };
  1573.  
  1574.  
  1575.   {
  1576.    double a = 1.5;
  1577.    double b = -1;
  1578.    double c;
  1579.    double s;
  1580.    double r_expected = 1.80277563773;
  1581.    double z_expected = -0.554700196225;
  1582.    double c_expected = 0.832050294338;
  1583.    double s_expected = -0.554700196225;
  1584.    cblas_drotg(&a, &b, &c, &s);
  1585.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 534)");
  1586.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 535)");
  1587.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 536)");
  1588.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 537)");
  1589.   };
  1590.  
  1591.  
  1592.   {
  1593.    double a = 1.5;
  1594.    double b = -0.1;
  1595.    double c;
  1596.    double s;
  1597.    double r_expected = 1.50332963784;
  1598.    double z_expected = -0.0665190105238;
  1599.    double c_expected = 0.997785157857;
  1600.    double s_expected = -0.0665190105238;
  1601.    cblas_drotg(&a, &b, &c, &s);
  1602.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 538)");
  1603.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 539)");
  1604.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 540)");
  1605.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 541)");
  1606.   };
  1607.  
  1608.  
  1609.   {
  1610.    double a = 1.5;
  1611.    double b = 0;
  1612.    double c;
  1613.    double s;
  1614.    double r_expected = 1.5;
  1615.    double z_expected = 0;
  1616.    double c_expected = 1;
  1617.    double s_expected = 0;
  1618.    cblas_drotg(&a, &b, &c, &s);
  1619.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 542)");
  1620.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 543)");
  1621.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 544)");
  1622.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 545)");
  1623.   };
  1624.  
  1625.  
  1626.   {
  1627.    double a = 1.5;
  1628.    double b = 0.1;
  1629.    double c;
  1630.    double s;
  1631.    double r_expected = 1.50332963784;
  1632.    double z_expected = 0.0665190105238;
  1633.    double c_expected = 0.997785157857;
  1634.    double s_expected = 0.0665190105238;
  1635.    cblas_drotg(&a, &b, &c, &s);
  1636.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 546)");
  1637.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 547)");
  1638.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 548)");
  1639.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 549)");
  1640.   };
  1641.  
  1642.  
  1643.   {
  1644.    double a = 1.5;
  1645.    double b = 1;
  1646.    double c;
  1647.    double s;
  1648.    double r_expected = 1.80277563773;
  1649.    double z_expected = 0.554700196225;
  1650.    double c_expected = 0.832050294338;
  1651.    double s_expected = 0.554700196225;
  1652.    cblas_drotg(&a, &b, &c, &s);
  1653.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 550)");
  1654.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 551)");
  1655.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 552)");
  1656.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 553)");
  1657.   };
  1658.  
  1659.  
  1660.   {
  1661.    double a = 1.5;
  1662.    double b = 1.5;
  1663.    double c;
  1664.    double s;
  1665.    double r_expected = 2.12132034356;
  1666.    double z_expected = 1.41421356237;
  1667.    double c_expected = 0.707106781187;
  1668.    double s_expected = 0.707106781187;
  1669.    cblas_drotg(&a, &b, &c, &s);
  1670.    gsl_test_rel(a, r_expected, dbleps, "drotg(case 554)");
  1671.    gsl_test_rel(b, z_expected, dbleps, "drotg(case 555)");
  1672.    gsl_test_rel(c, c_expected, dbleps, "drotg(case 556)");
  1673.    gsl_test_rel(s, s_expected, dbleps, "drotg(case 557)");
  1674.   };
  1675.  
  1676.  
  1677. }
  1678.